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ABSTRACT 

We present results of Swift optical, UV and X-ray observations of the afterglow of 
GRB 050801. The source is visible over the full optical, UV and X-ray energy range of 
the Swift UVOT and XRT instruments. Both optical and X-ray lightcurves exhibit a 
broad plateau {At/t ~ 1) during the first few hundred seconds after the 7-ray event. 
We investigate the multiwavelength spectral and timing properties of the afterglow, 
and we suggest that the behaviour at early times is compatible with an energy injection 
by a newly born magnetar with a period of a few tenths of a millisecond, which keeps 
the forward shock refreshed over this short interval by irradiation. Reverse shock 
emission is not observed. Its suppression might be due to GRB ejecta being permeated 
by high magnetic fields, as expected for outflows powered by a magnetar. Finally, the 
multiwavelength study allows a determination of the burst redshift, z = 1.56. 

Key words: Gamma-Ray Bursts. 



1 INTRODUCTION 

The Swift mission represents a major breakthrough for the 
study of Gamma-Ray Bursts (GRBs). Ne w bursts are de- 
tecte d by the Burst Alert Telescope (BAT, iBarthelmv et al] 
I2OO5I ). a coded-mask imager with higher sensitivity than that 
of the BATSE instrument onboard the Compton Gamma- 
Ray Observatory. Smft can begin observing a GRB with the 
X-ray telescope (XRT, Burrows c t al. 2005j) and Ult raviolet 
and Optical Telescope (UVOT. iRoming et al.ll2005h within 
~ 1 minute of the GRB onset, because of Swift's unprece- 
dented capability of fast and autonomous repointing. Swift 
has therefore enabled us to routinely explore these sources 
from the end of the prompt emission and the beginning of 
the afterglow, an epoch which was poorly known before this 
mission. The observations with the XRT and UVOT tele- 
scopes are addressing a series of key issues in GRB studies, 
such as the link between the 7-ray and the afterglow emis- 
sion, the mechanisms which produce them, the duration and 
the properties of the "central engine" , and the origin of the 
prompt optical emission. 

In this paper, we focus on GRB 050801. Among the bursts 
observed by Swift, this GRB was the first to have an early 



afterglow characterized by a peculiar phase of steady flux, 
both in the X-ray and in the Optical/UV, lasting a few 
hundred seconds. We will show that this distinguishing be- 
haviour might be explained if we assume that the burst 
ejecta receive an injection of energy from a magnetar, born 
in the event that caused the GRB, with an initial period of 
a fraction of milliseconds. 

GRB 050801c ould thus be ano ther example, in addition to 
GRB 051221 jFan fc Xull200^ b of a GRB with a magnetar 
as central engine. 



2 PROMPT EMISSION 

This burst was detected by the Swift satellite on August 1, 
2005 at 18:28:16 UT. The 7-ray lightcurve registered by BAT 
(15-350 keV) shows two peaks, with a duration tgo = 20±2 s 
(see Fig. [T|. Henceforth errors are at the 68% confidence 
level, unless otherwise stated. The spectral analysis was per- 
formed only in the 15-150 keV band, because this is the 
band in which BAT has the optimal sensitivity. We tried 
the Band, cut-olT powerlaw and simple powerlaw models to 
fit the data. None of them gave a fit significantly better 
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than the others. However, the parameters in the first two are 
poorly constrained. We have therefore adopted the results of 
the simple powerlaw model fit, which gives a spectral index 
P = —1.00 ± 0.16 with a reduced chi-squared of Xv = 1-43 
with 73 d.o.f.. The total fluence in the 15-150 keV band is 
2.88 X lO"'^ erg cm~^, which allows us to classify this burst 
as moderately faint. Interestingly, the second peak is softer 
than the first one, being visible only at energies lower than 
^ 50keV. In the following, we present the X-ray, optical and 
UV observations performed with Swift and we discuss the 
multiwavelength properties of the GRB 050801 's afterglow. 



3 SWIFT X-RAY AND OPTICAL 
OBSERVATIONS OF GRB 050801 

Following the BAT trigger, GRB 050801 was observed by 
the XRT and UVOT instruments. A bright, unknown X-ray 
and optical source was found inside the BAT error circle, 
at coordinates RA=13:36:35.4, Dec=-21:55:42.3 (as mea- 
sured with UVOT). This source subsequently faded, and 
it was therefore confirmed as the afterglow of this burst. 
Other prompt optical observations were successfully per- 
formed with ROTSE ( Rykoff et al 2006) , and further obser- 
vations were made bv lFvnbo et al.l l|2005ai lb[) . 6 hours and 30 
hours after th e burst onset respect ively. No radio afterglow 
was detected (jCameron et al ][200i). 

The UVOT began observing the field of GRB 050801 
52s after the BAT trigger. The first ~9.4s exposure, taken 
while the spacecraft was settling on the target, was taken 
through the V filter in photon-counting ("event") mode. 
Once the pointing had stabilised, a 100s V-band finding- 
chart exposure in combined image and event mode was 
taken. After this, the UVOT cycled through the colour fil- 
ters, firstly with 10s exposures, then 100s and finally 900s. 
The later exposures were taken in Image Mode (IM) only. 

A series of images were created from the settling and 
finding chart event lists. A single image was produced from 
the settling exposure, and the 100s finding-chart event list 
was used to produce 11 time-sliced images. The first 2 of 
these images have 5s exposure times, and the next 9 images 
are of 10s duration. Counts were extracted from the images 
in all filters using a 4" radius aperture at the position of the 
afterglow. A 25" radius region, offset from the source, was 
used to determine the background count rates. The source 
count rates were aperture-corrected to a 6" rad ius to ensur e 
compatibility with the photometry calibration (|Poolell2005h . 
and were corrected for detector dead-time and coincidence 
loss (i.e. lost counts due to multiple photons arriving in a 
single region of the detector within a single frame). The 
count rates were converted to magnitudes and fiuxes using, 
respectively, the zero-points and counts-to-flux conversions 
available in the Swift CALDB. The complete log of UVOT 
observations is shown in Table [T] 

Observations using the XRT began approximately 61s 
after the burst trigger. The first exposure, taken in the Im- 
age Mode (IM), did not show any source. Subsequently, the 
XRT observed in Windowed Timing (WT) mode for approx- 
imately 20s. This observational mode did not accumulate 
enough counts to produce any useful spectral information. 
The lightcurve was obtained from WT event data using Xs- 



elect with 53 arcsecond wide extraction slices for the source 
and background. 

The rest of the XRT data were obtained in Photon 
Counting (PC) mode. The initial PC exposure lasted ap- 
proximately 880s and was found to be piled up. Comparison 
of the observed point spread function of the source to the 
model XRT point spread function indicated that significant 
pile up is restricted to the inner 7"radius region. Therefore 
source spectra and lightcurves were extracted from an an- 
nulus with an inner radius of 7"and an outer radius of 71". 
Later data did not suffer from pile up and the extraction 
region was a circle of radius 70.7". The background extrac- 
tion region was taken as a circle of radius 189"for both the 
piled up and the non-piled up data. An additional fainter 
source was found within the extraction region and this was 
excluded for both the source and the piled up extraction re- 
gions with a radius of 2.5"radius. For the PC mode data we 
included events of grade 0-12, in which the charge is spread 
over ^ 4 pixels on the detector. The total exposure integra- 
tion time in PC mode was 2.3 x 10* s. For the first day of 
observations a hot column in the XRT was close to the cen- 
tre of the source, while at later times Smft was repointed 
to move the source away from this artefact. For spectral 
analysis we used the response matrix from the latest Swift 
cahbration database, CALDB 20051201, and the tasks XR- 
TEXPOMAP and XRTMKARF were used to produce effective 
area files, corrected for the hot column. The pileup correc- 
tion factor was also taken from the arf files. This correction 
factor was dependant upon the distance of the burst centre 
to the hot pixel columns. 

Figure [3 shows the X-ray lightcurve of GRB 050801 in 
the 0.2-10 keV energy band. The optical/UV lightcurve is 
shown in Fig.|3l the optical emission of GRB 050801 does not 
decline in the first ~ 200 s, after which it begins a system- 
atic decay. This behaviour is very similar for the lightcurves 
in all filters, indicating no variation in the optical spec- 
trum. We fitted the whole dataset of optical/UV datapoints, 
which have been normalized to the V band (Fig. |3l), with a 
broken powerlaw. In the following we use the convention 
F oc t~°'i/~^ , where F, t and u are the fiux, the time since 
the BAT trigger and the photon energy respectively. The 
best fit values of the parameters are: an initial decay slope 
Qo,i ~ 0.17±0.09, a break time to, 6 = 230± 15 s and a late 
decay slope 0:0,2 ~ 1-18 ± 0.05. The X-ray flux also lacks 
a strong decline (Alog_F < 0.5) in the first ~ 200 s, but 
decays after this with a rate similar to that of the optical 
lightcurve. By fitting this lightcurve with a broken power- 
law, we obtain an initial decay slope ax,i = —0.23 ± 0.24, 
a break time tx,b = 257^24 s, and a late decay slope of 
ax,2 = 1.22 ±0.04: all these values are consistent with those 
obtained by fitting the optical lightcurve, within 2a. In fact, 
we find that the optical-to- X-ray flux ratio remains remark- 
ably constant throughout the whole observation (see Fig- 
ure |4ll. In order to formally assess that, we deconvolved the 
intrinsic distribution of our data from the distribution of 
measurement errors. We followed the maximum likelihood 
method of Maccacaro et al. (1988) to obtain the best esti- 
mate of the intrinsic standard deviation aM of the optical- 
to-X ray flux ratio. We got aox,M = 0.02to oi, which is 
consistent, within 1.5 standard deviation, with 0. The 3a 
upper limit on intrinsic standard deviation is 0.06. 

We also note that the level of the optical fiux found by 
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ROTSE 22s after the trigger l|Rvkoff et al.ll2006l ) is consis- 
tent with the plateau we find 70-230s after the trigger, indi- 
cating that the optical flux was roughly constant for ~ 200s 
after the end of the prompt emission. 

The 0.2-10 keV spectrum obtained wfth XRT is well fit- 
ted with an absorbed powerlaw model; results are listed in 
Table [2] and Fig. There is no evidence for spectral evo- 
lution during the whole of the follow-up observation. This 
is confirmed by the fact that there is no change in the 
softness ratio, as shown in Fig. [6] By means of the max- 
imum likelihood method of Maccacaro, we found that the 
best estimate of the intrinsic dispersion of softness ratio is 
o"s,A/ = 0.075 ± 0.062, which is consistent with at I.Ict 
level. 

The similar behaviour of the X-ray and optical 
lightcurves, as well as the absence of any change in the spec- 
tral properties, suggests that these two bands both lie in the 
same spectral segment. In order to test this hypothesis, and 
to understand the spectral properties of the burst in more 
detail, we reconstructed the Spectral Energy Distribution 
(SED) between these bands, spanning ~4 decades in fre- 
quency. 

We performed a joint spectral fit of optical and X-ray 
data taken between 160 s and 970 s after the burst (i.e. 
where both X-ray and optical flux were highest). We flt- 
ted the spectra with a powerlaw model, which is widely as- 
sumed to be the spectral shape of GRB emission, and we 
have included the effects of photoelectric absorption using 
the phabs and zphabs models in XSPEC. In our flt, phabs 
is fixed at the Galactic value Nh = 7 x 10'^'' cm~'^ (Dickey 
& Lockman 1990), while zphabs is allowed to vary. The op- 
tical/UV datapoints were corrected for the Galactic dust 
extinction before reading into XSPEC, by using the Galac- 
tic extinction law reported in Seaton (1979). Such a multi- 
component model gives a rough reproduction of the whole 
spectrum, with = 38.4/27 (see Fig.[71). The best fit value 
of the optical to X-ray spectral index is l3ox ~ 0.85 ± 0.02, 
which is consistent with the spectral slope in the X-ray band, 
supporting a scenario in which the X-ray and optical data 
lie on the same spectral segment. 

Joint analysis of the optical and X-ray data has also 
allowed us to derive a photometric redshift z for the burst; no 
spectroscopic redshift is available. Extragalactic hydrogen 
will cause a dropout in the emission blueward of 91.2(1 + 
z) nm, which can be detected by Swift UVOT if z > 1. 
In the case of GRB050801, there is a dramatic difference 
in flux between the UVM2 and UVW2 fllters, (see figure 
7), and from the fit described above, we formally obtain 
z — 1.56 ± 0.06. We should, however, note that the spectral 
model we used includes the sharp cut-off at the Ly-limit, 
91.2(1 + z) nm, but does not include the Ly-Q forest, or Ly- 
Q absorption from the host galaxy of the GRB. Together, 
these two sources of Ly-Q absorption will erode the spectrum 
between 121.6 and 121.6(1 -I- z) nm, suggesting that the true 
redshift might be slightly below the value we determined 
and increasing the uncertainty on the redshift beyond the 
statistical errors on the model fit. At redshift ~ 1.6, Ly-Q 
from intervening systems is expected to absorb less than 10 
per cent of the broadband fiux (Madau 2005), and will only 
influence the photometry in the UV fllters of the UVOT. As 
this attenuation is smaller than the statistical uncertainty 
on the UV datapoints (15, 23 and 28 per cent for UVWl, 



UVM2 and UVW2 respectively) , we consider a more realistic 
uncertainty for the redshift to be 0.10. 

Despite these caveats, our photometric redshift should 
be close enough to the actual value to permit an adequate 
estimation of the energy emitted by GRB 050801 during 
its prompt emission phase. In calculating this energy we 
have taken the K-correction into account using the method 
of Bloom et al. (2001), extrapolating the flux to energies 
below the BAT threshold. Assuming that the energy index 
does not change for i5 < 15 keV, we flnd the energy released 
over the 15-150 keV band (and in the GRB restframe) to 
be 2.3 X 10^^ erg. This value puts GRB 050801 in the low 
energy tail of the GRB prompt energy distribution (O'Brien 
et al. 2006, Nousek et al. 2006). Assuming that the spectrum 
of the 7-ray emission does not change at higher energies, we 
calculate the total isotropic equivalent energy of this burst 
in the 1 - 10000 keV band to be 9.16 x 10^^ erg. 



4 DISCUSSION 

4.1 Onset of the afterglow 

4-1.1 Constraints on the energy injection mechanism. 

The main peculiarity of the afterglow of GRB 050801 is the 
broad plateau observed at all wavelengths from optical to 
X-ray during the early phase. As we can see from Figure [51 
the X-ray emission is practically constant for ~ 200 s, with 
a ratio between duration 5t and time t which is nearly 1. 
There is no evidence for sharp decays o r X-ray flares such as 
those observed in s everal Swift GRBs (|Burrows et al.ll2005l . 
iNousek et alll2006h . which are usually associated with "in- 
ternal shocks" between the ultra-relativistic shells emitted 
by the central engine. 

Among the various scenarios that can explain the ob- 
served plateau there are models in which the energy dis- 
tribution over the flreball surfa ce is non-uniform; the se in- 
clude the "patchy jet" models (iKumar fc Piran gOOCl ) and 
the "two component" jet models ( Zhang et al.ll2004^ . In the 
flrst case, the assumption is that the burst emission is re- 
leased in a jet with large fluctuations in angular direction, 
i.e. with a patchy distribution of Lorentz factor 7. Emission 
from different patches becomes observable at different times, 
according to their diverse values of 7"^ , and this can produce 
features which are observed in all bands almost simultane- 
ously (as in GRB 050801). It should however be noted that 
the presence of large fiuctuations of energy across the ejecta 
is still speculative, and that in order to give rise to early 
(< 300 s) bumps, a distribution over a very narrow angular 
scale (6 ~ 0.1°) is required. Moreover, theoretical predic- 
tions show that in this scenario the "bump" should follow a 
broad peak (see lZhang et al.ll2006al ). There is no indication 
of this in the available data. 

The "two component jet" model assumes the exis- 
tence of two components in the ejecta: one ultra-relativistic 
component powering the GRB, and a second -moderately 
relativistic- cocoon component. The cocoon is deceler- 
ated later than the ultra-relati vistic component because 
of its lower Lorentz factor (see Zhang fc Meszar"o3 |2004| . 
iMeszaros fc ReesI l200ll . iRamirez-Ruiz et all l2002f ). Its af- 
terglow emission would therefore start later and may give 
rise to a "bump" in the afterglow lightcurve. However, in 
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Filter 


T J 

^ rma 


T 

J- range 


T 


Mag 


V 


57 


52—61 


9.4 


1401 -1- n 07 


V 


67 


64—69 


5.0 


14.5 it 0.1 


V 


72 


69—74 


5.0 


14.8 ± 0.1 


V 


79 


74—84 


10.0 


1 K 1 o _l_ n nq 

±iJ.±'J _1_ VJ.L/i7 


V 


89 


84—94 


10.0 


1 K 07 -L n no 


V 


99 


94—104 


10.0 


1 K no + n OQ 


V 


109 


104—114 


10.0 


1 4 Qi -L n no 


V 


119 


114—124 


10.0 


1 c; riK -L n no 

±\j.\j'j _i_ 


V 


129 


124—134 


10.0 


1 4 8Q + n OS 


V 


139 


134—144 


10.0 


1 K 1 q + n OQ 


V 


149 


144—154 


10.0 


1 K 04 _l_ n nq 


V 


159 


154—164 


9.8 


1 4 02 + n 08 


V 


244 


2S9— 248 


9.8 


15.4 ± 0.1 


V 


328 


323-333 


9.8 


15.25 ± 0.09 


V 


412 


408-417 


9.8 


15.5 ± 0.1 


V 


497 


492 — 502 


9.8 


16.1 ± 0.1 


V 


581 


576 — 586 


9.8 


16.2 ± 0.1 


V 


666 


ggl_g71 


9.8 


16.8 ± 0.1 


V 


750 


745—755 


9.8 


16.6 ± 0.1 


V 


835 


SSO- 840 


9.8 


16.7 ± 0.1 


V 


919 


914—924 


9.8 


16.6 ± 0.1 


V 


5965 




899.8 


> 19.3 


V 


35314 


1 7583 — 53044 


2697.6 


> 20.6 


V 


1494g7 


1 08434-1 90500 


4211.8 


> 20.0 


B 


215 


210-220 


9.8 


15.32 ± 0.06 


B 


299 


294-304 


9.8 


15.57 ± 0.06 


B 


383 




9.8 


1 fi 1 Q -1- n 08 

-L"J. -LtJ _1_ VJ . WO 


B 


468 


463—473 


9.8 


1 6 29 + 08 


B 


553 


548 — 557 


9.8 


1 fi Q7 -1- n 08 

-Hj-iji _i_ w . wo 


B 


637 


R'?2 — 642 


9.8 


1 ft nq -1- n nq 


B 


721 


717—726 


9.7 


1 fi fin -1- n no 
J- "J. WW _i_ w . wo 


B 


806 


801—811 


9.8 


1 7 nf! -L n nq 

-L 1 .WW _1_ W.Wi7 


B 


890 


885—895 


9.8 


1 7 1 q -1- n nq 

i- 1 . ±iJ _1_ W . WC7 


B 


972 


Q70— 97'^ 


3.3 


17.1 ± 0.2 


B 


33617 


1 5767—51466 


3599.1 


> 21.1 


u 


201 


1 Qfi_205 


9.8 


1 4 QK -L n nfi 

-Lrt.iJW _1_ W . WW 


U 


285 


280-290 


9.8 


14.71 ± 0.07 


u 


369 


365-374 


9.7 


15.04 ± 0.08 


u 


454 


449-459 


9.8 


15.02 ± 0.08 


u 


538 


533-543 


9.8 


15.5 ± 0.1 


u 


623 


618-628 


9.8 


15.7 ± 0.1 


u 


707 


702-712 


9.8 


15.9 ± 0.1 


u 


792 


787-797 


9.8 


16.0 ± 0.1 


u 


876 


871-881 


9.8 


15.6 ± 0.1 


u 


961 


956-965 


9.8 


16.3 ± 0.1 


u 


12159 


11773-12545 


771.3 


19.30 ± 0.02 


u 


35301 


23342-47259 


2114.0 


20.99 ± 0.01 



Table 1. UV observations of GRB 050801 with UVOT in different time intervals after the BAT trigger. Errors are at 68% C.L. For those 
observations where no source was found with at least 2ct significance, we list upper limits at Scr. 



accordance with the theoretical models mentioned above, 
such a feature shou ld only be visible later than ~ lO'* s 
l|Zhang et aLlbOOeal l. 

The possibility we favour is that the plateau seen in 
both the X-ray and optical flux corresponds to the early 
afterglow, i.e. to the phase that follows the creation of an 
"external shock" after the fireball runs into the circumburst 
medium. The standard afterglow model predicts that, as a 
shell of ejecta interacts with the circumburst medium, it 
gives rise to a "forward shock", which propagates outward. 



In its simplest formulation, the model predicts the produc- 
tion of a broad and achromatic flare, rather than a plateau, 
both in the optical and X-ray bands. The non-decay of the 
afterglow emission of GRB 050801 before 250s after the trig- 
ger might then be explained by assuming that the central en- 
gine continues to inject radiative energy into the fireball for 
a few hundreds or even thousands of seco nds after the initial 
explo sion (late "central engine" activity. IZhang fc Meszaroj 
l2002h . As the emission from the central engine stops, the 
afterglow decay rate steepens up to more typical values. As- 
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Filter 


T J 

^ rma 




T 


Mag 


UVWl 


187 


182- 


-192 


9.8 


14.7 ± 0.1 


UVWl 


271 


267- 


-276 


9.8 


15.0 ± 0.1 


UVWl 


356 


351- 


-361 


9.8 


15.4 ± 0.2 


UVWl 


440 


435- 


-445 


9.8 


15.4 ± 0.2 


UVWl 


525 


520- 


-530 


9.8 


-LlJ . O _1_ "J. 


UVWl 


609 


604- 


-614 


9.8 


1 « n + n 2 


UVWl 


694 


689- 


-698 


9.8 




UVWl 


778 


773- 


-783 


9.8 


1 K 2 -1- n 2 


UVWl 


862 


858- 


-867 


9.8 


1 « q -1- n 

X.\J .iJ _1_ \J . \J 


UVWl 


947 


942- 


-952 


9.8 


1 f! 1 4 -1- n 2 


UVWl 


11316 


10866- 


-11766 


899.8 


1 Q 1 2 -1- n r)4 


UVWl 


34511 


22435- 


-46587 


2699.3 


> 20.6 


UVM2 


173 


168- 


-178 


9.8 


14.7 it 0.2 


UVM9 


257 


252- 


-262 


9.8 


14 9 it 0.2 


UVM2 


342 


337- 


-346 


9.8 


15.1 it 0.2 


UVM2 


426 


421- 


-431 


9.8 


15.9 it 0.3 




511 


506- 


-515 


9.8 


1 « 2 + 3 


UVM2 


595 


590- 


-600 


9.8 


16 + 03 


UVM2 


679 


674- 


-684 


9.8 


16 3 + 03 


UVM2 


764 


759- 


-769 


9.8 


16.5 ± 0.4 


UVM2 


848 


843- 


-853 


9.8 


17.8 ± 0.5 


UVM2 


933 


928- 


-938 


9.8 


17.4 ± 0.5 


UVM2 


6590 


6422- 


-6759 


337.1 


18.7 ± 0.1 


UVM2 


10415 


9966- 


-10864 


897.3 


1 q 4 -I- n 07 

J-c/.rt _i_ yj.yji 




39067 


21527- 


-56607 


2910.5 


> 20.6 


UVW2 


230 


225- 


-235 


9.8 


16.1 ± 0.2 


UVW2 


314 


309- 


-319 


9.8 


16.6 ± 0.3 


UVW2 


399 


394- 


-404 


9.8 


16.5 ± 0.3 


UVW2 


483 


478- 


-488 


9.8 


17.1 ± 0.3 


UVW2 


568 


563- 


-573 


9.8 


18.2 ± 0.4 


UVW2 


652 


647- 


-657 


9.8 


18.3 ± 0.5 


UVW2 


737 


732- 


-741 


9.8 


>17.9 


UVW2 


821 


816- 


-826 


9.8 


17.4 ± 0.4 


UVW2 


906 


901- 


-910 


9.8 


17.6 ± 0.4 


UVW2 


4524 


4188- 


-4859 


670.6 


20.2 ± 0.05 


UVW2 


34525 


16676- 


-52374 


3599.1 


>21.6 



Table 1. Continued. 



Segment 


Energy index /3 


Nh X 1022cm-2 


time interval (s) 


Plateau (piled up data) 


0.79 ±0.23 




90-300 


Whole piled up PC data 


0.70 ±0.12 


< 0.05 


90-975 


Non-Piled up PC data 


0.87 ±0.23 


< 0.06 


4097-52950 



Table 2. Best fit parameters for the X-ray spectrum of GRB 050801 in different time intervals after the BAT trigger. WT data analysis 
results are not shown because of the low statistics. Errors in /3 are at 68% C.L, while Njj uppe r limits are reported at 90 % C.L. and at 
^ = 0. The values of N^f represent the excess absorption above the Galactic value of 7 X 10^" ( ijPickev fc Lockman|[l99Ct )) 



suming that the luminosity of the central engine scales as 
L oc the predicted relation betwee n the decay and spec- 
tral slope is (jZhang fc Meszaros|[2002l '): 

a = {q/2 + l)l3 + q-l, (1) 

provided that the energy band we are dealing with is 
between the synchrotron peak frequency and the cooling 
frequency VclJ As we will show later, this is indeed the case 

^ Note that the signs of the indices have been reversed with re- 



fer both the optical and the X-ray observed frequencies. As 
for the decay index, we used the weighted average of ax,i 
and ao,i, which is a„ — 0.12 ± 0.09. For the spectral slope, 
we used Pox = 0.85 ± 0.02. From this choice we obtained 
g = 0.19 ±0.07. 



spect to the original paper bv lZhang fc MeszarosI ll2002h . in order 
to make them consistent with our convention (here both a and /3 
are positive). 
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This value of q is consistent with within 3 cr, and in- 
terestingly such a flat luminosity distribution is expected in 
models in which the gamma ray burst leads to the birth of 
a millisecond pulsar with an ultra-strong magnetic field, i.e. 
a so-called "magnetar" (see e.g. Zhang & Meszaros 2001, 
Dai & Lu 1998). If this is the case, several parameters of 
the newly-born compact object can be determined. We will 
discuss this intriguing scenario in more detail in § 14.1.21 
An observationally indistinguishable possibility is that the 
central engine activity is as brief as the prompt emission 
itself but, at the end of the prompt phase, the ejecta are 
released with different Lorentz factors (see Rees fc Me szaros 
1998, Sari & Meszaros 2000. IPanaitescu et al.ll2005h . The 
slowest shells would catch up with the fastest ones, once the 
latter have been decelerated by interaction with the circum- 
burst medium. The additional energy produced when the 
shells shock with each other would make the afterglow de- 
cay shallower than usual. In this case, assuming that the 
mass M of the ejecta follows the law M(> 7) (x 7", where 7 
is the shell's Lorentz factor, one can find an effective a value 
that mimics the effects of continuous energy injection with 
the luminosity law L ^ t'"^. Following iZhang et al., (|2006al ), 
this is: 



-(10-7g)/(2 + g) 



(2) 



and therefore the value q — inferred above is equiva- 
lent to s = —4.4 ±0.4, i.e. a very steep distribution. A steep 
distribution in the shells' Lorentz factors is required in order 
to have significant energy injection into the blast wave, with 
more energy carried by slow shells. Furthermore, the shell 
injection must occur over a suitably extended time interval 
in order to produce a flat plateau instead of a flare. 

As mentioned previously, both these interpretations 
hold if the band where observations are taken is positioned 
between Vm and vc- This can be inferred by means of the 
"closure relationships" , which link the value of the spectral 
index in a given band with the decay index of the flux in 
the same band, ac cording to the standard a fterglow scenario 
(for a review, see IZhang fc Meszarodliooi ). These relation- 
ships are different depending upon the the frequency range 
of the observed band, the density profile of the circumburst 
medium (constant as in the interstellar medium, or radi- 
ally decreasing outwards from the centre of the explosion as 
expected in the case of a wind from a massive GRB progeni- 
tor), and on whether the ejecta expand spherically or in the 
form of a jet. 

After the end of the plateau phase, in which (accord- 
ing to our interpretation) the behaviour of the afterglow is 
altered by late-time energy injection, one of these relation- 
ships should be satisfied if we assume that the afterglow 
behaves as theoretically predicted. In the case of the opti- 
cal and X-ray afterglow of GRB 050801, we find that the 
only closure relationship which is satisfied is a = (3/2)/3 
-which corresponds to the case of a spherical expansion of 
the ejecta into a constant density medium- and observing 
frequency v satisfying i>m < v < fc- This also allows us to 
obtain a value for p, the power-law index of the energy dis- 
tribution of the radiating electrons. On theoretical grounds, 
and for this spectral segment, it is a = (p— l)/2 which gives 
p = 2.70 ±0.05. 



4.1.2 Birth of a Magnetar? 

According to Zhang & Meszaros (2001), when the core of 
the stellar progenitor collapses into a pulsar, the luminosity 
emitted in the form of Poynting fiux, Lo, is given by 



10 BisP_.jRgergs 



(3) 



where we use the convention Qx = Q/10^. In the above 
expression, I, ilo, P, R and B are, respectively, the neu- 
tron star's moment of inertia, initial angular velocity, initial 
period, radius and initial magnetic field. The quantity Tem 
represents the time scale over which the emitted luminosity 
is roughly constant at the level Lo; following again the same 
authors, it is given by 



Tpm — 



_3£J_ 

B2R<^nl 



= 2.05 X 10^(/45-Br5 P-3i?^^)s. 



(4) 



After the time T^m, the model predicts a sharp cut off 
of luminosity L, which decreases as L t~ , i.e. q — 2. 
From Zhang fc Meszaros 2002, the afterglow dynamics can 
be affected only if q < 1. Therefore, after the Tem, we have 
no further modification of afterglow lightcurve due to energy 
injection. 

By multiplying Eq. |3] with |4l we can estimate the 
amount of energy that is produced by the newly-born neu- 
tron star and subsequently injected into the ejecta via dipo- 
lar spin-down. We get 



LoTem = 2.05 X 10"^ PI^ erg. 



(5) 



Our data also allow us to estimate the change in kinetic 
energy, Ek, during the "plateau" phase. In the case where 
the observing frequency is between the peak frequency and 
the cooling frequency, the kinetic energy released as a func- 
tion of the time is given by (Zhang et al., 2Q06b) 



Ek.; 



6.5 X 10-^1 



4/{p+3) 



d: 



,8/(p+3) 



X (l+^)-3{P-l)/(P+3)i3(p-l)/(p+3)^-4/{p+3) (g) 

y ^{-p-l)/(p+3) 4(l+p)/{p-|-3) -2/(p+3) 2{p-3)/{p-l-3) 
^ ^B,-2 E,-l "■ '^la 

where F{u = lO^^Hz) is the flux density at ly = 10^^ Hz, 
D is the luminosity distance, td the time in the observer's 
frame, fp is a function of the index p of the powerlaw en- 
ergy distribution of radiating electrons, en is the fraction 
of the energy that is contained in magnetic fields, es is the 
fraction of the energy associated with radiating electrons, n 
is the particle density in the surrounding medium and f is 
the frequency at which we are observing. Again, subindices 
indicate normalized quantities {Qx = Q/10^), while td is in 
units of days. 

For a GRB at z=1.56, assuming the cosmological pa- 
rameters Ho ~ 70, n — 0.3 and A — 0.7, we find 
_D = 4 X 10^* cm. For the following parameters we tak e 
fiducial values of e_B = 0.3 ([Freedman fc Waxmann|[200ll ). 
es = 0.03, n = 10 (|Berger et all |2003| ). u = 10^^ Hz, 
p = 2.7 and fp = 0.1 (the function fp is plotted in Zhang 
et al. 2006b). Inserting the time duration of the plateau 
[Aid = (250 - 70) s « 2.1 x 10"^ d] into Eq. [gl and using 
the 0.2 - lOkeV flux, F ^ 4x 10"^°erg cm"^s"\ we can 
infer the increase in kinetic energy during this phase. This 
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gives /S.Ek ~ 1.4 X 10^"^ erg. Finally, under the simple as- 
sumption that all the electromagnetic energy emitted by the 
pulsar has been completely converted into the kinetic energy 
of the ejecta (we discuss more realistic scenarios below), we 
have Ai?K = LoTem and Eq. Ogives the neutron star period, 
Po « 0.4 X 10"^ s. 

We can now use the value of the spin period to 
estimate the emitted luminosity and the pulsar's magnetic 
field strength. In the cosmological rest frame of the GRB, 
the plateau duration is (250 - 70) x (1 + z)^^ « 70 s, 
which gives a good approximation for the quantity Tem. 
From Eq. [S] we therefore obtain a luminosity Lq ~ 2 x 10^^ 
erg sec~^ and, from Eq. O a stellar magnetic field of 
i? ~ 1.4 X 10^^ G (using a neutron star radius of 10 km). 
This value is in the range expected for ultramagnetized 
neutron stars ("magnetars"). 

The existence of a plateau in the lightcurve also requires 
the pulsar to have a significant effect on the energetics of 
the afterglow. This means that the energy injected by the 
pulsar before the beginning of the afterglow, which we set 
at a certain time To, must be a sizeable fraction of the total 
kinetic energy that is left for the ejecta after the prompt 
emission phase, E'l^: 

LqTq > E'k . (7) 

To test whether this condition is met, let us consider the 
energy radiated in 7-rays during the prompt emission. 
The radiative efficiency of a GRB, defined by the relation 
T)^ = E-y/[E^ + E'l^), gives a measure of how efficiently the 
GRB dissipates the total energy into radiati on during the 
7-ray emission phase. On theoretical grounds (|Guetta et al.l 
[20o3) this efficiency is not expected to be very high, with 
typical values of ~ 0.2 which are in agreement with pre- 
Swift results (De Pasquale et al. 2006) . More recent analysis 
indicates much higher efficiencies, up to 0.7 — 0.8, although 
these values are hard to accommodate within the standard 
"internal shock" emission mechanism of GRB prompt emis- 
sion. Assuming a conservative value of 77^. = 0.2, and taking 
E^ = 9.16 X 10^^ erg for GRB 050801 (see § [3|, we ob- 
tain E'lf = 3.6 X 10^^ erg. This means that the condition 
LoTo ^ E'li is met for To > 20s. 

Unfortunately, a precise estimate of To, the time of the 
beginning of the afterglow, is not possible. Based on our 
data, we can only state that To must be less than ~ 70 s, 
when Swift began its follow-up observations. Comparing 
with the BAT lightcurve, we suggest that the afterglow be- 
ginning might be identified as the time of the second peak in 
the prompt emission, which is ~ 5 s after the first peak. As 
mentioned above, this peak is broad and clearly softer than 
the first peak, which is to be expected if the first peak were 
due to internal shocks and the second to afterglow emission. 
Such behaviour has already been observed in the prompt 
emission of several GRBs, su ch as GRB 970228 (for a re- 
view see iFrontera et"aLl I2OO0I ). It is worth noting that, if 
the afterglow begins about ~ 5 s after the trigger, the op- 
tical emission detected by ROTSE can be associated with 
the afterglow and connected with the flat optical lightcurve 
recorded by UVOT later on (as observed). This comes with 
the caveat, though, that since we do not know the behaviour 
of the GRB in X-rays in the interval 20- 70s after the trigger. 



we cannot be sure that the emission in this period is actually 
due to the forward shock. 

We can, however, discuss the possibility that the energy 
injection, which shapes the early afterglow curve, begins at 
the time of first ROTSE observation, i.e. 20s after the trig- 
ger. In such a case, we obtain, from Eq. 5, that the energy 
injected is 1.26 higher than previously calculated. This in 
turn decreases the required initial period P by a factor of 
1.12, i.e. we obtain P — 0.36ms. From Eq. 3, we can infer 
that the magnetic field would lower by a factor 1.26, so we 
have B ~ 10^^ G. As we can see, there is not a big change 
in the values of the parameters derived. 

Finally, we shall briefly discuss the consequences of a 
slightly incorrect estimate of the redshift, as previously pre- 
sented. If the true redshift is slightly above the value we 
used, the value of Lq and T^m estimated previously should 
be increased and diminished respectively. For small correc- 
tions of z, the product of these two parameters (see Eq. 5) 
does not vary largely. In addition, P has a weak dependence 
on this product. Therefore, we do not expect a large varia- 
tion of P. For example, if the true redshift were 10% lower 
than the value we reported, the value of P should increase 
only by 3%. Conversely, if the true redshift were 10% higher 
than that we found, P should be 3% lower. Similarly, we 
do not predict a noticeable variation of the magnetic fleld. 
From Eq. 3, we can infer that the value of P15 would change 
by ~ 1% for a 10% variation of the redshift. In the same 
way, the condition expressed by Eq. 7 is still satisfled in 
our model, for small changes in z without any considerable 
change of To- 

4.1.3 Critical issues 

The discussion presented in § 14.1.21 is based on a number 
of assumptions regarding the parameters which describe the 
physical properties of the GRB. We note, however, that the 
dependence of the kinetic energy of the ejecta on es, es, 
and n is not very strong, and so different values are not ex- 
pected to signiflcantly affect our conclusions. On the other 
hand, the assumption that all of the EM energy radiated by 
the pulsar is converted into the kinetic energy of the ejecta 
is more questionable. In a more realistic scenario, we would 
expect that only a fraction of this EM energy is converted 
in this way; in order for our scenario to be consistent, the 
efficiency cannot be lower than ~ 30%. We certainly do not 
expect the conversion efficiency to be too low, however, since 
the relatively low-energy emission from the magnetar is un- 
able to penetrate the Shockwave. The Poynting fiux energy 
of the pulsar can therefore only propagate outwards by con- 
version into the kinetic energy of the ejecta. The constraints 
regarding the conversion of the EM energy may be some- 
what further eased if we assume that the ejecta outflow is 
moderately beamed. While such a condition would not be 
in contradiction with the observations or theory, it would 
further reduce the energy contribution required from the 
pulsar. 

Other uncertainties relate to the way in which the pulsar 
energy injection is estimated. We assume that the magnetic 
fleld is constant during the phase in which energy injection 
is associated with the rapid pulsar spin down, while in a 
more realistic scenario the coupling between rotational and 
magnetic evolution should be accounted for. Moreover, the 
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emitted EM energy is calculated under the simple assump- 
tion of dipolar emission in vacuum, which is obviously a 
crude approximation during the first phases of neutron star 
formation and evolution. Dynamo effects, the presence of 
magnetospheric currents or multipolar components of the 
magnetic field may complicate this simple picture. 

Our overall conclusion is that, to a first order approxi- 
mation, EM energy injection from a rapidly spinning mag- 
netar (a neutron star rotating with a period of a few tenths 
of a millisecond and possessing an ultra-strong magnetic 
field of order ~ fO^^ G) can plausibly explain the observed 
plateau in the X-ray and optical lightcurves of GRB 050801. 
Although we have not directly derived estimates of the 
burst's physical parameters (eB,e_B, density of surrounding 
medium) , the assumption of canonical values for these quan- 
tities is an adequate approximation for the purposes of our 
calculations. 

It is inst ructive to c o mpar e GRB 050801 with the short 
GRB 051221 iFan fc Xul (|2006l '). In this latter event, a flat- 
tening occurred in the X-ray lightcurve between 3 x 10^ s 
and 2 x 10* s after the trigger. Optical data on this burst 
are rather sparse, and so cannot help constrain the afterglow 
physics. In the context of short GRBs, which are thought to 
arise from the coalescence of a binary system composed of 
two compact objects (such as two neutron stars or a neu- 
tron star and a black hole) , it is difficult to imagine that the 
fall-back accretion of part of the material onto the central 
compact remnant can, hours after the coalescence, pump 
the emitted energy up to 10^^ — 10^^ erg. Another scenario 
was therefore put forward in which two neutron stars co- 
alesce and form a magnetar. In this hypothesis, the flat- 
tening of the lightcurve was interpreted as the signature of 
the underlying magnetar, whose energy injection can only 
significantly influence the afterglow lightcurve after several 
kiloseconds post-trigger. According to these authors, a mag- 
netar with an initial period of 10~^s and magnetic field of 
10^* G can explain the observed behaviour. In the case of 
GRB 050801, we infer values for the period and magnetic 
field which are, respectively, shorter and higher. This is be- 
cause, for this burst, the magnetar energy injection must oc- 
cur over a shorter time interval, whilst the luminosity must 
be higher. It should be noted, nevertheless, that the differ- 
ences in these parameters between the two GRBs are within 
~ 1 order of magnitude. 

4.2 Absence of a reverse shock? 

In the standard afterglow model, it is expected that the 
formation of a forward shock is accompanied by that of a 
less energet ic "reverse shock'' whic h moves inward through 
the ejecta (|Meszaros fc Ree3ll997l ). While the emission of 
the forward wave is predicted to peak in the X-ray band 
(at least at the early stages of the afterglow that we are 
considering here), the emission of the reverse shock should 
peak in the IR-optical because of the higher density of the 
ejecta that it crosses. 

In the case of GRB 050801, we do not observe any opti- 
cal ffares, and the optical emission always follows the X-ray 
emission, with the ratio of optical to X-ray ffux remarkably 
constant (see Fig. This suggests that the reverse shock 
emission is suppressed. A possible reason may be a very high 
magnetic field in the ejecta, in which case most of the energy 



is carried by the field itself and it is not converted in t o rad i- 
ation in the shells. According to lZhang fc Kobavashil l|2005h . 
for bursts with short duration (i.e. fgo ,$20 s), suppression 
of the reverse shock takes place if the parameter a, defined 

as 

De,2 

""^-A J' (8) 

where B'^ and n are the magnetic field and the den- 
sity of the ejecta respectively, and rrip is the proton mass, 
is larger than 100. The a parameter can be interpreted as 
the ratio between the energy contained in the magnetic field 
and that in the baryonic outflow. Since the value of B'^ is 
strongly dependent on the assumptions made about the po- 
sition of the shock radius and the local particle density, it 
cannot be robustly determined. If internal shocks occur at 
small radii, even a relatively low magnetic field may pre- 
vent the reverse shock from forming. We can, however, give 
an order-of-magnitude estimate of assuming an ejecta 
mass of 10~* — 10~^ M©, a shock radius at the onset of the 
afterglow of r^h ~ lO^^cm and a spherical shell of depth 
c&t with &t ~ tgo ~ 20 s. This gives a particle density of 
n ~ 10** cm^'^ in the shell, which in turn implies a magnetic 
fleld of B'^ ~ 10* G. Another observational consequence of 
a strongly magnetized reverse shock is the broadening of 
its emission profile (Zhang fc Kobayashi 2005). This phe- 
nomenon, together with the suppression of the peak emis- 
sion, would make the reverse shock emission even less evi- 
dent. It should be noted that the hypothesis of a high mag- 
netic field in the ejecta may be consistent with the powering 
of the wind by a magnetar. Thus, the non-detection of a re- 
verse shock component in the emission of this burst may be 
in agreement with the identification of the central engine of 
GRB 050801 as a magnetar. 

Other scenarios are, however, also plausible. For exam- 
ple, we cannot exclude the possibility that, for unusually 
high density ejecta, the reverse shock emission would have 
peaked in the far-IR rather than in the IR or optical. In such 
a case, it would have made a negligible contribution to the 
emission detected by UVOT, which is consistent with our 
observations. 



5 CONCLUSIONS 

Swijt GRB 050801 shows remarkable spectral and temporal 
features, in both the X-ray and optical/UV bands. Instead 
of rapidly decaying fiux within 100-200s after the trigger, as 
seen in the majority of bursts, the GRB 050801 emission in 
the X-ray and in the optical has a plateau from ~ 70 s up 
to ~ 250 s after the BAT trigger, followed by a more normal 
decay slope. It is possible that the plateau actually extends 
back to ~ 20 s after the trigger, which could explain why the 
first prompt optical measurement, performed by ROTSE-III 
at this epoch, is consistent with the UVOT data. 

We find that the relatively fiat X-ray lightcurve might 
be caused by a late energy injection by the "central engine" 
to the expanding shell, in the form of either Poynting fiux 
or the shock dissipation that occurs when late shells catch 
up with earlier ones. In the case of energy injection due to 
Poynting fiux, we have shown that dipolar spinning down 
emission from a newly born magnetar of initial period Po — 
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0.4 s and magnetic field B ~ 1.4 x 10^^ G can account 
for the X-ray flux observed during the plateau phase and 
for its duration. Other models, which involve uneven energy 
distribution in the ejecta or internal shock emission are likely 
to be ruled out by the fact that the plateau appears at early 
time and its emission is remarkably smooth. 

The energy injection model requires that the afterglow 
has already started by the time of the Smft measurements. 
This hypothesis is supported by the fact that the optical-to- 
X-ray flux ratio during the plateau is consistent with that 
observed at late times, when the afterglow has certainly be- 
gun. We suggest that the afterglow onset might actually be 
associated with the second peak detected by BAT. This peak 
is very soft, as expected if it is connected with the afterglow 
emission, and it occurred only ~ 5 s after the trigger. 

Theory predicts that the forward shock in the circum- 
burst medium, responsible for the afterglow emission, should 
be accompanied by a reverse shock that moves inwards 
through the ejecta. The reverse shock emission should be 
brief and peaked in the optical. On the other hand, the con- 
stancy of the optical to X-ray flux ratio observed during 
this event suggests that no reverse shock emission has taken 
place, at least during the Swift observations. The suppres- 
sion of the reverse shock emission might suggest the presence 
of highly magnetized ejecta, in which most of the energy is 
carried by the magnetic field rather than the shock. 

Finally, a joint fit of the X-ray and optical data has 
allowed us to confirm that data in both bands lay on the 
same spectral segment, between the synchrotron peak and 
the cooling frequency, and to determine the spectral index 
with high precision. Analysis of the late-time (> 300 s) after- 
glow emission, based on the application of closure relations, 
suggests that the GRB fireball had a spherical expansion in 
a constant density environment. More importantly, we have 
been able to determine the redshift of the burst as z = 1.56. 
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Figure 1. The GRB 050801 BAT lightcurve. From top to bottom; lightcurves in the 15-25 keV, 25-50 keV, 50-100 keV, 100-350 keV, 
and 15-350 keV energy bands. Units on the Y axis are counts per second per detector. 
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Figure 2. The GRB 050801 X-ray afterglow observed with Swift XRT in the 0.2-10 keV energy band. 
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Figure 3. Swift UVOT count rate of GRB 050801 in all filters, renormalized to the V band. We also show the ROTSE R band points, 
renormalized with the same criterion. 
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Figure 4. Ratio of the Swift XRT and UVOT count rate (renormalized). The ratio is defined as (Cx-Co)/(Cjf +Co)i where Cx and 
Co are the count rates in the XRT and a given UVOT filter respectively. 
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Figure 5. The GRB 050801 X-ray spectrum during the piled-up phase (time interval 90-975s, top points) and the non piled-up phase 
(since 975 s, bottom points). 
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Figure 6. Softness ratio of X-ray data, defined as (Co. 2-2 — C2-io)/(Co.2-2 + C2-10), where Co. 2-2 and C2-10 are the count rate in 
the 0.2-10 IjeV band and in the 2-10 keV band, respectively. 
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Figure 7. Spectral Energy Distribution of GRB 050801, iromSwift XRT and UVOT data. The best-fit model is also plotted. In this 
figure, the model and XRT datapoints were obtained by unfolding the best fit model in XSPEC. The optical/UV fiux densities were 
obtained by multiplying the UVOT countrates by the standard countrate-to-flux conversion factors (Poole 2005). 



